Superfluid density and penetration depth in Fe-pnictides 
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We consider the superfluid density p s (T) in a two-band superconductor with sign-changing ex- 
tended s-wave symmetry (s + ) in the presence of non-magnetic impurities and apply the results to 
Fe-pnictides. We show that the behavior of the superfluid density is essentially the same as in an 
ordinary s-wave superconductor with magnetic impurities. We show that, for moderate to strong 
inter-band impurity scattering, ps(T) behaves as a power-law T" with n ~ 1.6 -r 2 over a wide 
range of T. We argue that the power-law behavior is consistent with recent experiments on the 
penetration depth X(T) in doped BaFe2As2, but disagree quantitatively with the data on LaFePO. 



Introduction. Recent discovery of iron-based pnictide 
superconductors instigated massive theoretical and ex- 
perimental research effort aimed at unveiling fundamen- 
tal properties of these materials. Both oxygen containing 
'llll's materials (La,Nd,Pr,Sm)FeAsO, and oxygen free 
'122's (Ca,Ba,Sr)Fe2As2 have high potential for the ap- 
plications and may create a breakthrough in the field of 
superconductivity (SC)pi 

One of the central and still unsettled issues is the sym- 
metry of the SC gap. An ordinary s-wave superconduc- 
tivity due to phonons has been deemed unlikely because 
of too small electron-phonon coupling^ suggesting that 
the SC pairing is of electronic origin. Electronic struc- 
ture of pnictides shows pairs of small hole and electron 
pockets centered at (0,0) and (%, 7r), respectively in the 
folded Brillouin zonei 2 -^ Most of parent compounds dis- 
play an SDW order with momentum at, or near (7r,7r), 
and an early scenario was the pairing mediated by an- 
tiferromagnetic spin fluctuations^ For pnictide geome- 
try, this mechanism yields an extended s-wave gap which 
changes sign between hole and electron pockets but re- 
mains approximately uniform along either of them (an 
s + gap). The s + gap symmetry has been found in weak 
coupling studies of two-band^ and five-band^ "g-ology" 
models of interacting low-energy fermions. Some other 
studies, however, found a gap with extended s-wave sym- 
metry in the unfolded Brillouin zone£ Such gap has no 
nodes on the hole Fermi surface (FS) , but has four nodes 
on the electron FS, like a d-wave gap in the cuprates. 
The uncertainty arises from the fact that in pnictides 
there is a competition between the inter-pocket interac- 
tion with large momentum transfer and intra-pocket re- 
pulsion. When inter-pocket interaction is stronger, the 
system likely develops a sign-changing s— wave gap with- 
out nodes; when intra-pocket repulsion is stronger, the 
system develops an extended s— wave gap with nodes to 
minimize the effect of intra-pocket repulsion— 

^From experimental perspective, the situation is 
also unclear. Andreev spectroscopy^ and ARPES 
measurements^ are consistent with the gap without 
nodes. NMR and Knight shift measurements^ were orig- 
inally interpreted as evidence of a gap with nodes, but 
it turns out that the data can be fitted equally well by 
a dirty s + superconductor^^ The situation is further 
complicated by the fact that the two hole FS are of dif- 



ferent sizes, and have different gapsJ^ 

To truly distinguish between s + gap and a gap with 
the nodes one should go to very low temperatures. Re- 
cently, two groups reported measurements of the pene- 
tration depth A(T) down to 10 -2 T c . The data seem to 
point into different directions. Ames group reported the 
data on Co-and K-dopcd BaFc2As2 (Ref. fl5T) and demon- 
strated that down to the lowest T and for all dopings 
the T dependence of A(T) = A(0) + AA(T) can be fitted 
by AA(T) oc T 2 . In SmFcAsOi_ x F 2 i& and PrFeAsOil 
penetration depth seems to have an exponential temper- 
ature dependence at low T, consistent with the gap with- 
out nodes. At the same time, in another 1111 material, 
LaFePO, the Bristol group has found AA oc T 12 down to 
the lowest temperatures^ FS in this material has been 
reconstructed from magneto-oscillation measurements^ 
and consists of weakly corrugated small-size cylinders 
making it unlikely that either hole or electron FSs extend 
to the points where s + order parameter changes sign. 

The penetration depth in the clean limit was consid- 
ered in Ref. for an s + gap and in Ref. [2l] for several 
other gap symmetries. For a gap without nodes, AA(T) 
is obviously exponential at low T, for a gap with nodes 
it is linear in T. In this paper we discuss to what ex- 
tent the existing data for A(T) can be described by a 
dirty SC with an s + gap symmetry, and non-magnetic 
impurities. We argue that the T 2 behavior observed by 
Ames group in 122 materials can be fitted over a wide 
temperature range for various dopings and in this respect 
the results for the penetration depth are not in conflict 
with ARPES and other measurements which show SC 
gap without nodes. The existence of the two different 
gaps on the two hole FS makes the agreement between 
experiments even better. On the other hand, the linear 
in T behavior observed in LaFePO is not reproduced, and 
it is more likely that the gap in this material has nodes, 
as suggested in Ref. H. 

For a conventional s-wave SC the effect of non- 
magnetic impurities on the temperature dependence of 
A(T) is small and mostly irrelevant for all T. For s + su- 
perconductors, the situation is qualitatively different be- 
cause inter-band impurity scattering mixes hole and 
electron states with opposite values of the order param- 
eter ±A and in this respect should be pairbrcaking and 
act in the same way as a magnetic impurity in a conven- 
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tional s-wave superconductor. Consequently, scattering 
by non-magnetic impurities in s + SC affects T c , the den- 
sity of states, and the temperature dependence of the 
penetration depth. Over some range of IV/A, the be- 
havior at the lowest T is still exponential, however when 
Ttt/A becomes larger than a critical value, superconduc- 
tivity becomes gapless, and the exponential behavior dis- 
appears. 

Method. The London penetration depth A(T) scales 
as 1/y p s (T), where p s {T) is the supcrfluid density. The 
latter is, up to a factor, the zero frequency value of the 
current-current correlation function and can be written 
in the forrn^ 
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where p s o is the superfluid density at T = in the ab- 
sence of impurities. The integrand in Eq. ([T]) is defined 
in terms of impurity- renormalized Matsubara energy, <D m , 
and the superconducting vertex A rn . In an s + supercon- 
ductor the order parameters on the hole (c) and electron 
(/) FS pockets are related, A^ = — A£, = A m and in 
Born approximation 

iQ m = iu> m - T g c (uj m , A TO ) - T n g f (Cj m , A TO ) , (2a) 
A m = A + r / c (£ m ,A m ) + r^(£ m ,A m ), (2b) 

where u> m = irT(2m + 1), To = ittiiN f\u$\ 2 and = 
imiNp^^ 2 are the intra- and inter-band impurity scat- 
tering rates, respectively {u lT are impurity scattering 
amplitudes with correspondingly small, or close to 7r = 
(tt, 7r), momentum transfer), A is the SC order parame- 
ter, and functions g c ^ and f c ^ are £— integrated normal 
and anomalous Green's functions for holes and electrons: 



g c = g f = 
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Since the /-function has opposite signs in two bands, 
r,r has the same effect on anomalous self-energy as the 
scattering on magnetic impurities in an ordinary s-wave 
superconductor. Following the customary path one may 
introduce r) m = uj m /uj m and A m = A m /r] m that satisfy 



Vm = 1 + (r + T n ) ■ 
A m = A(T) - 2T* 
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The order parameter A(T) is determined by the self- 
consistency equation 



A(T) = V sc nTj2f c (Q m ,A m ) = ttT^ 



V sc A r 



u m V m 1 m 

(5) 



where V sc is the s + coupling constant and A is the ul- 
traviolet cutoff. Notice that the last expression contains 
A rn and bare Matsubara frequencies oj m . Eqs. d2J)— (J5j) 
can be extended to the case when the gaps on hole and 
electron FS have different magnitudes. 

Solutions of the system of Eqs. (|4b[) and ([5]) give the 
values of A(T) and A m . In particular, Eq. (|4b[) is an 
algebraic equation (valid at any T) which expresses A m 
in terms of A. The latter itself depends on rV, because 
the self-consistency equation (J5J contains A m . Without 
inter-band scattering (r^ = 0) we have A m = A = Ao = 
1.76 T c o, where T c o and Ao are the BCS transition tem- 
perature and the T = gap in a clean superconductor. 
For 7^ 0, A m differs from A, and A differs from Ao. 
Converted to real frequencies, Eqs. (|4b[) and (O yield a 
complex function A(oj). For 2r T > A, A(ui = 0) van- 
ishes, i.e., superconductivity becomes gapless<22 At the 
critical point 2T 7T = A, A(uj) cx (— iw) 2//3 at small u), at 
larger T n , A(ui) — —i const ui + 0{io 2 ). 

Results. We express the results using dimension- 
less parameter £ = T n /2TrT C Q. For equal gap magni- 
tudes and 2T 7r /A < 1, y = A/A is the solution of y = 
cxp[— 7re 7 C/y], where 7 « 0.577 is the Euler constant^ 
At a given T a gapless superconductivity emerges, when y 
becomes smaller than 4£e 7 , i.e., for ( > (1/4) exp[— (7 + 
7r/4)] w 0.064. The transition temperature obey a 23 i 24 
ln(T c /T c0 ) = # (1/2) - * (1/2 + 2(T c0 /T c ), where V(x) 
is the di-Gamma function, Fig. QJi. T c decreases with 
C and vanishes at Q cr = e" 7 /8 w 0.07 (T^/Ao = 1/4). 
For 0.064 < ( < ( cr , A(T, lo) cx ioj for small u>, includ- 
ing T = 0, and thus even T = zero-energy density 
of states becomes finite. (At the onset, at £ = 0.064, 
T c w 0.22 T c0 , and A(0) = 0.46 A ). The ratio 2A(0)/T C 
increases with £ and reaches 7.2 at the onset of the gap- 
less behavior and 8.88 at ( = ( cr (Ref. [24]). A large value 
of 2A(0)/T C is often attributed to strong coupling^ but, 
as we see, can also be due to impurities. 

In terms of auxiliary A m and rj m , 
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In general, the value of p s (T = 0) and the functional 
form of p s (T) depend on both ITV and Tq because To is 
explicitly present in Eq. ([6]) via t] m given by Eq. ((4a]). 
Impurity scattering amplitude is a decreasing function of 
momentum transfer, and, in general, To ^> T v . Since we 
are interested in ~ A, we have T S> A and 



A; 



A: 



(7) 



where B = ■np s0 /{T + T n ). We see that r only affects 
the overall factor B, and all non-trivial T dependence 
comes from frequency and temperature dependence of 
A m . 

Several results for p s (T) given by (J7J) can be obtained 
analytically. First, near T c , p s (T) cx A 2 (T) cx T c - T, 
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FIG. 1: (color online) (a) Suppression of T c by inter-band 
scattering in a two-band (A, —A) model; (b) superfluid stiff- 
ness p s (T) in a dirty s + superconductor for fixed intra-band 
impurity scattering To/2irT c o = 3, and various inter-band 
scatterings £ = T^/2tvT c o; and (c) low-T plot of p s vs. T 2 
showing near n — 2 power-law around onset of gapless regime. 

i.e., 

where p s (T = 0) is the actual zero-temperature value 
of p s . In a clean BCS superconductor B = 2. In the 
present dirty case (T ^> T c , T^) B(Q is non-monotonic 
in C and equals B(( ->■ 0) « 2.65, B(C = 0.064) w 
1.67, S(C « Ccr) = 2.03. This implies that a lin- 
ear extrapolation of p s from T « T c to T = still 
yields a significantly larger value than the actual p s (0). 
Second, at £ < 0.064, the T dependence of p s {T) re- 
mains exponential at low T, p s (T) oc e~ A ' u=0 '' T with 
A(w = 0) = A [l - (C/Ccr) 2/3 ] 3/2 ,_but at the onset 
of gapless superconductivity, when A(w) oc (— iw) 2 / 3 , 
we have p s (T) oc T 5 / 3 . Finally, in the gapless regime 
0.064 < C < Ccr, we found p s (T) oc T 2 at low T. 

To obtain p s {T) at arbitrary T, we self-consistently 
solved the gap equation (J5]) together with impurities ([2]) 
and Green's functions ([3]) numerically, found A(T) and 
A m , substituted them into Eq. |T} and obtained p s (T). 
We present the results in Fig. [1] for several values of £. 

We see that, once the inter-band impurity scattering 
increases, the range of exponential behavior of p s (T) pro- 
gressively shrinks to smaller T, Fig. [TJj. Outside this low 
T-range, the temperature dependence of p s strongly re- 
sembles T 2 behavior, see Fig. QJc). The T 5 / 3 behavior at 
the onset of gapless superconductivity is hard to see nu- 
merically, as this power is confined to very low T, while 
for slightly larger T the behavior is again close to T 2 . 
Overall, the behavior of the superfluid density in a rela- 
tively wide range of £ is a power-law T n with n reasonably 
close to 2 down to quite low T. At the same time, we 
didn't find conditions under which the superfluid density 




FIG. 2: (color online) The fits to experimental data for 
BaFe 2 As 2 (RefQJ). We use only low-T data as at higher T 
the experimental A(T) may be influenced by sample geometry 
and fluctuations, (a) The data for electron-doped BaFe2As2 
for optimally doped (x = 7.4%) and overdoped (x = 10%) 
samples can be fitted reasonably well using (A, — A)-model; 
(b) The fit of data for two Co-doped (a; = 7.3%, 10%) and 
one K-doped samples with a phenomenological extension of 
the presented two-band model to the case of four gaps. We set 
gaps to be (A, -A; A, — Aj,) with A h = A/3, A/2. In this 
case the pairbreaking parameter does not need to be large. 
The fitting values of T = penetration length are large but for 
x — 7.4% sample we show a fit with inclusion of Fermi-liquid 
effects that reduces this parameter to A F l (0) = Ao ~ 190 nm, 
in agreement with experimental values. 



would be linear at low T. 

Comparison with the data. Judging by the value of 
2A(0)/T C fi£ the material with the least amount of inter- 
band impurity scattering is SmFeAsOi-^F^, where T c ~ 
55A'. In this compound it is difficult to expect a large 
£, since the observed exponential BCS-likc behavior of 
p s {T) at small T (Ref. \lw is consistent with extended 
s-wave gap and weak inter-band impurity scattering^ 

The data for electron- and hole-doped BaFe2As2 
(Ref. [HI) is fitted in Fig.[5^,b. The measured p s (T) scales 
approximately as T 2 , which is similar to behavior shown 
in Fig. [ljc). Left panel is the fit assuming that the gaps 
on two electron FSs and two hole FSs are A, — A; right 
panel is a more realistic fit in which we assumed, guided 
by ARPES data,ii that the gaps on the inner hole and 
the two electron FS are the same, but the gap on the 
other FS is 2-3 times smaller. 

The values £ = 0.04 — 0.06 used in these fits corre- 
spond to T c /T c q ~ 0.6 — 0.3 which is consistent with 
the values of T c ~ 10 — 30A' in this material, if we as- 
sume that T c o in the clean case is roughly the same as 
in SmFcAsO. The curves shown in Figs. [2^,,b represent 
the best fits, but we emphasize that we do not need to 
adjust £ to get a T 2 behavior - it persists over a range 
of C, see Fig fib. However we note that A(0) used in the 
fits is rather large compared with the experimentally ob- 
tained values ~ 200 — 300 nm^ We suggest that this 
discrepancy may be due to the omission of Fermi-liquid 
effects. The qualitative argument, supported by numcri- 
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cal estimates, is as follows. Assume, by analogy with the 
cuprates, that fermion-fermion interactions renormalizc 
lu — > luZ^, where Z u is a decaying function of frequency, 
and further assume that Z^ ks 1 at energies compara- 
ble to A so that it does not affect the relation between 
A(0) and T c . Then we find that the low-temperature de- 
pendence f(T/T c ) of the penetration length is rescalcd 
AA/Ao ~ f(ZoT/T c ) and the value of the fitting parame- 
ter Ao = A FL (0) decreases compared to what is obtained 
without Fermi liquid effects. This Z factor is particularly 
relevant to heavily underdoped regime where it increases 
because of SDW fluctuations. We believe this is the rea- 
son why we have to use very large A(0) = 2800nm to fit 
the data. At larger dopings, Z is smaller, but according 
to ARPESj" Z ~ 2 in optimally doped Bai_ a; K x Fe2As2. 
For / ~ T 2 , from the analytic reasoning we get the ef- 
fective Ao four times smaller than A(0), reducing it to 
Ao ~ 150 — 400 nm, in the range of what is experimentally 
extracted. The numerical analysis confirms this and for 
Z = 2 we show a fit for one of the electron-doped samples, 
that gives a reasonable value for the zero-temperature 
penetration length. We therefore conclude that the pen- 
etration depth data for 122 material can be fitted by a 
model of a dirty s + superconductor. Note in this regard 
that the data£& show that T c is almost insensitive to the 
value of residual resistivity. This was interpreted as the 
argument for a conventional s-wave gap. We note that 
this is also the case for extended s-wave gap as the dom- 
inant impurity scattering is intra-band scattering, con- 
trolled by To, which affects residual resistivity but does 
not affect T c . 

We attempted to fit the data for LaFePO using this 
approach but the fit fails for all £ and one has to assume 
unrealistically large A(0) to get even mediocre agreement 



with the data. From this perspective, it is likely that the 
linear in T behavior of A(T) in LaFePO is not related to 
impurities but rather is the consequence of the fact that 
the gap in this material has nodes. 

Conclusions. In this paper we considered superfluid 
density p s (T) in a multi-band superconductor with sign- 
changing s + - symmetry in the presence of non-magnetic 
impurities and applied the results to Fe-pnictides. We 
showed that the behavior of the superfluid density is es- 
sentially the same as in an ordinary s-wave superconduc- 
tor with magnetic impurities. For a moderate inter-band 
impurity scattering, p s (T) over a wide range of T behaves 
roughly as T 2 and crosses over to exponential behavior 
only at very low T. When superconductivity becomes 
gapless at T = 0, the T 2 behavior extends to the lowest 
T. We argue that this power-law behavior is consistent 
with recent experiments on penetration depth A(T) in 
hole and electron doped BaFe2As2, but we find that the 
present model does not explain the data for LaFePO. 

Several modifications of the model may improve the 
comparison with the experimental data. First, an ex- 
tension beyond Born approximation may further flatten 
the T dependence of A(T), although recent study for a 
different model 28 didn't find much changes beyond Born 
approximation. Second, we assumed that the surface of 
a superconductor sample is homogeneous. Defects on the 
surface may also modify the temperature dependence of 
the penetration depth. 
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